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ABSTRACT 

Analysis was performed to estimate the statistical 
visibility time spans of earth orbiting satellites as seen 
simultaneously by a groiond station and a ship. The analysis 
covers topics such as time average population, average population 
times and also the percentage visibility times for a given 
number of satellites. These results are useful for specific 
communications satellite applications. Numerical results are 
obtained for various configurations of ground station and ship. 
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SECTION 1 


INTRODUCTION 


This report is concerned with the analysis of the number and also 
the time of satellites mutually observed by both a ground station and a 

ship. Unlike the relatively simple case of a single observation sta- 
tion for which the region of observation is the volume bounded by a 
cone, the present more'' complicated case has a region of observation de- 
termined by the intersection of two cones. This region has a volume 
determined only by the separation distance between the ground station 
and the ship; but it also has a directional property determined by the 
relative position of the ship with respect to the ground station. 

Because the analysis becomes extremely complex, it is necessary to make 
certain simplifying assumptions. 

The first assumption is that the satellites presently orbiting the 
earth may be broadly classified into a few categories. This simplifica- 
tion is supported by the fact that^^^ since 1977 approximately 635 sat- 
ellites have been launched and these may be characterized as in Table 1.1 




Table 1.1 



Class 

J^verage 

Period 

Average 

Inclination 

Average 

Altitude 

Number 

I 

100 min. 

o 

O 

00 

800 km 

440 

II 

12 hr. 

60° 

20,000 km 

106 

III 

24 hr. 

0° 

36,000 km 

57 

IV 

Others 



32 

Thus , 

instead of having 

to deal with the volume 

of the region 

of observe' 


tion, the analysis deals with the areas at the various altitudes. In 
this analysis, only Class I and II satellites are considered. Class III 
satellites are considered separately because they are geosynchronous. 
Class IV satellites are irregular and will not be considered at all. 


(1) NASA, Satellite Situation Report, Volume 21, Number 1, 
February 28, 1981. 
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The second assumption is that within each of the two categories con- 
sidered, the satellites have circular ^orbits which are uniformly distrib- 
uted in terms of equatorial crossing and, moreover, the satellites are 
also uniformly distributed along the orbital arcs* 

Section 2 deals with the derivation of the number density of sate- 
llites in this statistical distribution. Section 3 deals with the deter- 
mination of the common region of observation of both the ground station 
and the ship. Section 4 is concerned with the computation of the time 
average population of satellites within the mutual region of observation. 
Section 5 briefly discusses the computation of the average population 
times of these satellites in the same region. Section 6 summarizes the 
results of this sfudy for Class I and II satellites. 

Readers who are strictly interested in the numerical results may 
go directly to Section 6 and omit the intervening sections which deal 
with the mathematical analysis. 
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SECTION 2 


STATISTICAL DESCRIPTION OF ORBITING SATELLITES 


2.1 Distribution Function 

Consider a statistical description of a system of N satellites as 
previously described in which the circular orbits are uniformly distributed 
in terms of equatorial crossing and the satellites are uniformly dis- 
tributed along the orbital arc. Consider Figure 2 . 1 which illustrates 
a given orbit with inclination i. Let 9 be the latitude, cj) be the right 
ascension measured from the equatorial crossing, and a be the orbital arc 
measured also from the equatorial crossing. 



<p 


Equator 


Figure 


2.1 


Consider Figure 2.2 which illustrates the area element dA at the equator 



Figure 


2.2 



It is obvious that dA^ is given by 


dA = r d<l) d0 
o o o 


( 2 . 1 ) 


where r is the radius of the orbit. Let f denote the density of the satel- 

o 

lites at the equator. Then, the number dN of satellites contained in dA^ 
is given by 


dN = f dA 
o o 


( 2 . 2 ) 


As these satellites move to latitude 0 and right ascension the corre- 
sponding area dA is then given by 


dA = r COS0 d(f>d0 


(2.3) 


and the density f is then obtained from 


dN = f dA 


(2.4) 


Substitution of Equations (2.1) - (2.3). into (2.4) yields 


£ d(p ■ d0 
0 0 o 


f = 


COS0 d^d0 


(2.5) 


However, from Figure 2.1, we obtain the following spherical trigonometric 
formula 


sin0 = sin i sin a (2.6) 


so that at latitude 6 we have 

cos0 d0 = sin i cosa do (2.7) 


6^5 



and at the equator we have 


d0 = sin i do 
o o 


Moreover, it is easily verified that we also have 

d<}> = d<},^ 

da = da 


( 2 . 8 ) 


(2.9) 

( 2 . 10 ) 


Substitution of Equations (2.7) - (2.10) into (2.5) yields 

f 


f = 


( 2 . 11 ) 


cos o 


which states that the density is inversely proportional to the cosine of 
the arc length. 

Next, we obtain the equatorial density as follows: 


- 0 max 

2. 

^ 2 

= 


f r 

^ 0 min 

J o 


f 

, 0max 


= 2..^ j 

0 min 

f cos0d0 


/2 


= 2TTr^ 

-7T/2 

f sin i 

= 27Tr^ 

r 7t/2 

f sin i 


0 


-7T/2 


o 2 2 

f sin 
o 

i 

= 2tt r 


( 2 . 12 ) 


in which Equations (2.7) and (2.11) have been used. 
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Substitution of Equation (2.12) into (2.11) yields 

N 


f 


0 2 2 
2n r 


sin i cos a 


(2.13) 


which expresses the density in terms of the total number, the radius, the 
inclination and the orbital arc. However, it is more convenient to obtain 
an expression in terms of latitude than orbital arc. This is accomplished 
as follows: Using the identity 


2 

cos a 


1 - sin a 


(2.14) 


and also Equation (2.6), we obtain 


”2 2 

sin i cos a = /(sin i - sin 0) 


(2.15) 


so that Equation (2.13) becomes 

N . 

f= ^2.16) 

o 2 2 77 . 2. , 

/7T r /(sin 1 - Sin 0) 

2.2 Angular Separation 

The system of N satellites under discussion is considered to be uni- 
formly distributed in terms of equatorial crossing and also along the or* 
bital arc. It is easily verified that the angular separations betx^een the 
satellites are given by 

(2.17) 

^ (2.18) 
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Figure 3.1 
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The geocentric conical angle a may be obtained as follows: For the triangle 

OGH, we have the sine formula 


(r + h) sin a 
e 

sin (tt - B) 


(3.1) 


and the cosine formula 


= r ^ + (r + h)^ - 2r (r + h) cos ot 
e e e e 


(3.2) 


Substitution of Equation (3.1) into (3.2) and use of the identity 


. 2 , 2 
sxn a = 1 - cos a 


(3.3) 


yield the following quadratic equation for cos a 


2 / r 

cos a - 2/ e 


2 

sin 6 cosa + 


Ir + h 
e 


r \ .2 2^ 

e \ sin 6 - cos 3 


^r + h 
e 


= 0 (3.4) 


whose solution is given by 


cosa = 


sin 3 + cos 


B / 


+ hJ 


1 - 


. 

Sin 3 


r + hj 
e 


(3.5) 


It may be verified that the physically acceptable solution is the one which 
yields the smaller angle a, i.e., the one with the positive sign in Equation 
(3.5). The other solution yields the larger angle a which results in a cone 

going into the earth, which is thus rejected. 

3.2 Boundaries of Intersection Region 

Consider Figure 3.2 which illustrates a ground station G and a ship 
S, and also the region of observation common to both of them. Let C be 
the central point of the great circular arc GS, and y the angle GOG. 
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Then, from spherical trignometry, the arc length 2y between G and S is 
given by 

cos 2 Y = cos 6 g cos ( 3 . 6 ) 

and the inclination K of S with respect to G is given by 


63 = sin K sin 2y 


(3.7) 


Next, consider Figure 3.4 which illustrates the boundaries R and L 
of the intersection region Q. It is to be noted that these boundaries 
are not arcs of great circles, but are arcs of small circles. 



In order to obtain expressions for the boundaries R and L, it is conven- 
ient to consider the arc GS as the equator in an oblique coordinate 
system. First, consider the curve R. 
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Let 0' be the latitude and \p" be the longitude of a point with respect 
to G. Then, from spherical trigonometry, the equation of the curve R is 
given by 

cosa = COS0* cost);" (3.8) 

However, if denotes the longitude measured from C, then we have 

ip” = ilj' + r (3.9) 

and Equation (3.8) becomes 

cosa = cos 0' cos (ij^' + y) (3.10) 

which is the equation for the boundary R in the oblique geographical system 
having C as the origin of latitude and longitude. Similarly, the equation 
for the boundary L is given by 


cosa = cos O' cos (ip' - y) 


(3.11) 


The points of intersection of the curves R and L are given by 
(f' = 0, 0' = e'max) and (^' =0, 0' = 0'min) where 




max 



( 3 . 12 ) 


0 ' . = - 0 ’ 

min max 


(3.13) 
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3.3 Regular to Oblique Geographic Transformation 


Let /r = (x,y,z) denote the coordinates of a point in the regular 
geographic system, and (x* , y’ , z*) denote the corresponding coordinates 

of the same point in the oblique geographic system. Figure 3.5 illustrates 
the angular rotations to accomplish the necessary coordinate transformation. 





(3.15) 


6^13 


where the symbols s and c respectively denote the sine and cosine func 
tions of the argument which appears as the subscript. Next, it is also 
noted that r and r' may be respectively expressed in terms of their 
latitude and longitude as follows: 


X = r c.c, 
0 


y = ^"4- 


z = r s. 




rs 


0 ' 


(3.16) 


(3.17) 


Thus, Equations (3.14) - (3.17) may now be used to express the oblique 
latitude and longitude in terms of the regular ones. The final results are 
given by 


s c^s , + c 
K 0 l/) K 0 


(3.18) 


tan (ip' + y) 


C c^s , + s, s . 
K 6 Ip K 6 


^9^tp 


(3.19) 


3.4 Oblique to Regular Geographic Transformation 

Next, to obtain the regular latitude and longitude in terms of the 
oblique ones, we proceed, as follows: We note that 

r = A^r’ (3.20) 

T 

where A denotes the transpose matrix of A. 



Then, proceeding as before but now using Equations. (3.20) and (3.15) - 
(3.17), we obtain 








+ c 


S c„, 
K e 




+ c 


^ 0 * 


(3.21) 


tan ip = 


+ Y) 


- s 


' 0 ' 


‘'e' ''(</'’ + Y) 


(3.22) 
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SECTION 4 - 


TIME AVERAGE POPULATION 


4.1 Exact Folmulation 


Let N denote the number of satellites (time average population) 
within the common domain of observation 0. It is obviously given by 



f dO 


(4.1) 


where the density f is given by Equation (2.16) and the element of area 
dn is given by 

dfi = r^ COS0 difi d6 (4.2) 


It appears that the above integral may be trivially expressed in terms of 
the regular latitude 0 and longitude i|/ as follows: 


N COS0 di|) d0 ■ (4.3) 

2 2 2 
2ir^ (sin i - sin 8) 


where ip (9) denotes the expression obtained by solving for \p in terms of 9 
R 

using the equation for the R curve given by Equations (3.10), (3.18), and 
(3.19), and similarly for ij;^(9) in terms of the L curve. Not only is this 
process difficult, but it is noted that the integral on the RHS of Equa- 
tion (4.3) may not even be valid or, worse yet, amenable to numerical eval 
uation even in principle. This point becomes evident by combining Figures 
3.3 and 3.4. It is possible that the location of the ship S with respect 
to the ground station G can give rise to the case where, in performing the 


0 max 




“a ■ 


9 rain 
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integration with respect to ij/, the process does not take place from the L 
curve to the R curve and, furthermore, in performing the integration with 
respect to 0, the process also does not take place f rom 0 min to 0 max . 

This difficulty can be circumvented by writing the element of area dJ^ 
as follows 

2 

dfi = r COS0' dij^' d0' (4.4) 


so that the Integral becomes 


"q = 


/ 


0 ^max 


0 *min 


(0*) 

r K 




N COS0 ^ d\|) * d0 * 


27t (sin i - sin 0) 




■J / 


cos0 ' d\lj ' d6 ' 


0'min ip' (e') 

li 






'•] 


(4.5) 


in which Equation (3.21) has baen used. It is to be noted the integra** 
tion will always proceed from the L curve to the R curve, and the 9* 
integration will always proceed from 6* min to 9^ max. 

4 . 2 Approximate Formulation 

An approximate formulation may be obtained by going back to the original 
Equation (4.1) which may be used to yield the following 


N.; = f ^ 
u ave 


(4.6) 
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where 


■# 


Q = (if) r" COS0' di|^' d0' 

Q 


2 /''J'r (o') 

= Ar^ / ^ COS0’ dif.' d.0*^ 

Jo Jo 


0'max 

4r'" / COS0' 

0 


/ 


cos 

1 cosa \ 

•" 

- Y 


Icose'i 



d0' (4.7) 


in which Equations (3.10) and (3.12) have been used. This integral may be 
evaluated numerically once the relative position of the ship S is speci- 
fied. The average value of f to be used may be obtained by averaging the 
4 values at the mid-points on the axes of symmetry of 0. These, in turn, 
may be obtained by averaging the values at the center C and those at the 

extremities P illustrated in Figure 3.4. Thus, we may write 
1 


f 

ave 


1 

8 


f(P^) + f(P2) + f(P3) + f(P^) + 4f (G) 


(4.8) 
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SECTION 5 - AVERAGE POPULATION TIMES 

5*1 Exact Formulation 

First consider Figure 2.1, for which we may write the following spher 
ical trigonometric formulas 

sin 0 
cos o 

where the relevant quantities have already been previously defined in 
Section 2. Next, consider Figure 5.1 which illustrates the ground station 
at G, when the satellite crosses the equator at N^. Subsequently, when 
the satellite has moved to latitude 0, the rotation of the earth has taken 
the ground station to the point G. 



sin 1 sin a 
cos 0 cos d) 


(5.1) 

(5.2) 


Figure 5.1 

Then, it is obvious that the following relation holds for both direct 
(i < 7 t/ 2) and retrogade (i > tt/ 2) orbits 


where 


u -f 4) = -f ujO 


U) ^ — 

p 

e 


(5.3) 

(5*4) 
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y “ longitude of satellite crossing measured from ground station 

(\> = right ascension of satellite measured from equatorial crossing 

if; = longitude of satellite measured from ground station 

a = orbital arc of satellite measured from equatorial crossing 

(jL) = ratio of satellite orbital period P to earth rotational period P 

e 

Substitution of Equation (5.3) into (5.2) yields 

COSO = COS0 cos (if; + wo - y) 

Equations (5.1) and (5.5) express the latitude and longitude 
of the orbital arc. Symbolically, we may write 

0 = 0 (a ; i) 

If; = If; (o; i, y) 

In turn, these equations may be substituted into Equations (3.18) and 
(3.19) to yield expressions for the oblique latitude 0’ and longitude 
if;* in terms of orbital arc a. Thus, we have 

0' = 6' (6, >p; <, y) = 0* (o; i. y) (5.8) 

= 4/' (6, 4^; <. y) = 4»' (o; i, M, K, y) (5.9) 

which constitute 2 equations in the 3 unknowns 9', ip' and a. If we wish 

to determine the point of intersection with the R curve bounding one side 
of the common region of observation Q, we also have Equation (3.10) which 
is 

cos 'u = cos 0' cos ('<p' + y) (5.10) 

Substitution of Equations (5.8) and (5.9) into (5.10) yields a com- 
plicated equation for a which may then be solved numerically to obtain 


(5.5) 
in terms 

(5.6) 

(5.7) 
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the value •? = corresponding to the intersection point. Next, to 

obtain the point of intersection with the L curve, we have Equation 
(3.11) which is really Equation (5.10) with Y replaced by -Y . Thus, 
the same process may be used to obtain the value 0=0 corresponding 
to the intersection point. Thus, the population time T of the satellite 
within the region is exactly given by 


Let denote the value of P which corresponds to the orbit passing 
through the central point C. The above process is first performed with 
a value P = P^ + A<|)^ where A(J)^ is a random number in the range 0< AcJ)^ < A<j) 
where A(|) is given by Equation (2.17) which is ° 


(5.12) 


The process is then repeated with values (p + nAc|)) where n = ±1, ±2,.. 
until no more orbits intersect the regionii . After this, the entire 
above process is then repeated with other random values of Ac)) . The 

average population times are then obtained by averaging the re^sults of 
all these processes. 


5.2 Approximate Formulation 

Consider Figure 5.2 which illustrates the spherical triangle formed 
by the equator, the meridian and the arc length of the central point C 

measured from the ground station G. This spherical triangle is fixed on 
the rotating earth. 





Figure 5.2 
6-21 



Then, we 


have the following spherical trigonometric formulas 


sin0 

= sinK 

siny 

(5.13) 

COSY 

= COS0 

cost|^ 

(5.1A) 

siny 

sinX = 

sini(; 

(5.15) 


which may be used to compute the latitude and longitude of C and also 
the angle X the arc GC makes with the meridian through C. 

Next, consider Figure 5.3 which illustrates the spherical triangle 
formed by the equator, the meridian and the orbital arc of a satellite 
just passing through the point C. This spherical triangle is fixed on the 
celestial sphere, which is inertial. 


c 





Figure 5.3 

Then, we have the following spherical trigonometric formulas 


sinS 

= sin i sino 

(5.16) 

GOSCr 

= COS0 cosc 

(5.17) 

sine 

sine = sine 

(5.18) 


which may be used to compute the orbital arc, the right ascension and 
also the angle C the arc NC makes with the meridian through C. 
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However, because of the earth’s rotation, the satellite’s ground 
track does not really make an angle ^ with the meridian through C. 
Rather, it is deflected through an angle 5 which is, in general, 
given by 

tan 5 = h) cose cosC (5.19) 

1 - o)cos i 

where co is defined by Equation (5.4). (It may be verified that this de- 
flection causes direct orbits to be more inclined and retrograde orbits 
to be less inclined as viewed by their ground tracks.) Thus, the angle 
between the satellite ground track and the meridian at point C is given 
by ( 5 -^), as shown in Figure 5.3. Next, consider Figure 5.4 which illu- 
strates the inclination r| of the orbital arc with the oblique equator 



Figure 5.4 


Hence, it is seen that we have 


n = A - C + ^ for ascending orbits 


n==A 4 -^- 5 - 7 T for descending orbits 


(5.20) 


It may also be verified that these equations are algebraically valid for 
both direct and retrograde orbits. 
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Next, consider Figure 5.5 which illustrates the case of a satellite 
just passing through the point D whi^:h is displaced by Ay from the central 
point C. This corresponds to an orbit whose equator crossing is displaced 
by A(() from the point N. 


D 



Figure 5.5 

Then, using spherical trigonometric formulas, it may be shown that Ay 
is related to A(j) by the following equation 


tan (y + Ay) = 


tan i sin (ip - 4) -f A(j)) 


tan i cos < cos - cj) 4* A<p) - sin k 


( 5 . 21 ) 


Thus, by replacing y with (y+Ay ), Equations (5.13) -■ (5.20) may be used 

to compute the inclination n of the new orbital arc with the oblique equa- 
tor. It may be verified that Equation (5.21) is also algebraically valid 
for both the cases of i > K and i < K. Moreover , it is also valid for 
both direct and retrograde orbits. Furthermore, it is valid for arbitrary 
finite differences A(|) and Ay , but considerable care must be exercised 
when taking the inverse tangent to obtain (y + Ay) in the correct quadrant 
corresponding to the increment A(|) . 



Up to this point, no approximations have been made. It is now assumed 
that the satellite's ground track is an arc of a great circle lying within 
the region and making an angle X] with the oblique equator GS. Figure 
5.6 illustrates the cases of orbital arcs passing through the points C and D. 



Figure 5.6 

Now, it is possible to write the following two approximate relations for 
the orbit passing through the central point C 

sin 6' “ sinn sin a* 

cos 0* = cos 0' cos 

where O' is the arc length measured from the oblique equatorial crossing. 
These two equations are the crude analogs of Equations (5.8) and (5.9) of 
the exact case. If we wish to determine the point of intersection with 
the R curve, we also have Equation (3.10) which is 


(5.22) 

(5.23) 


cosa = cos 6' cos (i/^* +7). 


(5. 24) 
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However, instead of substituting Equations (5.22) and (5.23) into (5.24) 
to yield a complicated equation for a’, it turns out to be the case 
that an algebraic equation can be obtained involving sin 0*. This is 
accomplished as follows: From Equations (5.22) and (5.23), the following 

auxiliary equation is obtained 


sin i/;* = 


tan e' 
tan n 


(5.25) 


Equation (5.24) is then written as 

cosa = cos 0* (cos cosy - sin 


= cos a’ cosy - cos 0' tan 0 ' 

tan n 



sin 0 ^ 
tan n 


sin y) 


siny 


siny 


or equivalently 

2 2 

cos y/^sin n - sin 0*) = sinn cosa + cosn siny sin d\ (5.26) 


By squaring both sides of this equation, it is obvious that a quadratic 
equation is obtained involving sin 0* . After much simplification, it may 
be shown that we have 


sin 0 * 
sin n 

A little consideration will reveal that for the intersection point with 
the R curve, it is necessary to retain only the positive sign in the above 
equation. Thus, this expression corresponds to the value at 0 “ 9 r* 


2 2 2 

- cosa siny cosn ± cosy (sin a - sin y sin n ) (5,27) 
(1 - sin Y sin n ) 
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However, from Equation (5.22), it is seen that the value is 
given by 

sine = Sinn sin (5.28) 


Consequently, we have 


sin o ’ “ ~ cosa siny cosn cosy /p 


sin a 


.2 . 2 , 

sin Y sin n) 


2 2 

(1 - sin Y sin n) 


(5.29) 


Next, to obtain the intersection point with the L curve, we have 
Equation (3;11) which is really Equation (5.24) with y replaced by -y . 
Thus, the same process may be used to obtain which can be shown to 
be given by retaining the negative sign in Equation (5.27). Consequently, 
we have the following result 

sin a * 

JLi 


2 2 2 
= cosa siny cosn -cosy/(sin a - sin y sin^ 


nl 


(5.30) 


2 2 

(1 - sin y sin n) 


which states that o ^ ~ expected (only for the case of the 

Lt K 

orbit passing through the central point C) . Thus, the population time 
T of the satellite within the region Q is approximately given by 


T 


2i\ 



(5.31) 


which is the crude analog of Equation (5.11). 

Next, to obtain the intersection point between the R curve and 
the orbit passing through the point D, a little consideration will re- 
veal that it suffices to replace y by (y + Ay) and also use the 
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corresponding value of n and then repeat the process above for computing 
a ' as given by Equation (5.29). However, to obtain the intersection 
point between the L curve and the orbit passing through D, a little more 
caution is now necessary. It now suffices to replace y by (-y + Ay) and 
also use the corresponding value of q and then repeat the process above 
for computing a„' but now retain the negative sign. This result yielding 
the value of is no, longer trivially the negative of a^' as for the 
Special case of C. 

The above process is first performed with a random value in the 
range 0 £ < A(}> where A(j) is given by Equation (2.17) which is 

A(() - ^ (5-32) 

/n 

The process is then repeated with values (Acf)^ + nA(|)) where n = +1, +2,... 
until no more orbits intersect the region Q* After this, the entire above 
process is then repeated with other random values of A({)^. The average 
population times are then obtained by averaging the results of all these 
processes. 

Finally, it must be mentioned that in order to insure that the cor- 
rect segment of the R circle (see Figure 5.6) is identified to yield the 
desired intersection point as given by the general analog of Equation 

(5.29) , a little consideration will reveal that we must have n in the 

range -tt/2 < n 7r/2. Thus, if n is outside this range, we must accordingly 
add to or subtract tt from n. Similarly, the same procedure applies to 
insure the identification of the correct segment of the L Circle to yield 
the desired intersection point as given by the general analog of Equation 

(5.30) . Furthermore, considerable thought will reveal that this assign- 
ment of the n range not only correctly gives the desired intersection 
points for orbits crossing the oblique equator within the observation re- 
gion n, but also for the case of equatorial crossings outside it for a 
range of A'v exceeding tt/2 measured from the central point C. The rea- 
sons for this are not apparent and, at first sight, it would seem that 
this assignment of n values outside the region Q leads to incorrect an- 
swers. But this is not so because of the manner in which the inverse 
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trigonometric functions are assigned their principal values. Thus, 
Equations (5.29) and (5.30) contain many subtle featpres in logic which 
automatically combine to yield, in mutual accord, the correct intersection 
points regardless of the equatorial crossing. In particular, additional 
consideration will reveal that it is only necessary to consider equatorial 
crossings such that the orbits intersect the oblique meridian through the 
central point C at an oblique latitude 9' not greater than 6* given by 



This corresponds to a range Ay* given by 



(5.33) 


(5.34) 


(5.35) 


or equivalently 


Ay* = min ( ti/ 2, sin~^ /~ < 


2 2 
cos y - cos g 


(5.36) 


2 2 2 2 
cos y - cos a + cos a tan n 


It is not difficult to see that if an orbit intersects the oblique 
equator outside the range Ay* and also eventually intersects the observa- 
tion region Q, then this orbit would already have been counted as lying 
within the acceptable range on the other side of the central point. 
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SECTION 6 


RESULTS FOR ORBITING SATELLITES 


6-1 Average Population Time Computations 

Computations were performed, except for minor modifications, according 
to the method discussed in Section 5.2 to obtain the average population 
times for Class I and II satellites. The representative values of para- 
meters used are shown in Table 6.1. 

Table 6 ♦ 1 


Quantity 

Class I 

Class II 

Period P (minutes) 

100.9 

717.9 

Inclination i (degrees) 

74.0 

63.9 

Altitude h (km) 

800 

20,178.5 

Number N 

400 

100 


The value of 6, the conical observation angle at the earth’s surface, is 
taken to be 80° for both the ground station and the ship. The ground sta- 
tion is taken to be at the origin of longitude and latitude while the ship 
is taken to be at various values of longitude and lafitude 0^ only in the 
first quadrant. It may be verified that for locations of the ship in the 
other quadrants, the corresponding results may be obtained by taking mirror 
reflections about the primary axes. 

After the average population times t have been obtained, the results 
were divided by the characteristic time T defined by 

( 6 . 1 ) 

to obtain the number of satellites visible to both the ground station and 
the ship. (T is the time for a satellite to travel the intra-satellite 
distance Ao where Ac is given by Equation (2. 18)..- The relevant results 
for Class I and II satellites are respectively summarized in Figures 6.1 
and 6.2, each of which was obtained by averaging the results using ^ given 
by Equation (5.T9) and those using ^ = 0. 
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Note: (1) Numbers in the boxes denote 

the number of satellites vi- 
sible for the percent of time 
indicated. 

(2) Numbers below the boxes de- 
note the relative longitude 
and the absolute latitude of 
the ship. 
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Figure -6.1 - Results for 100 Minute Orbiting Satellites 
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Note: (1) Numbers in the boxes denote 

the number of satellites vi- 
sible for the percent of time 
indicated. 

(2) Numbers below the boxes denote 
the relative longitude and the 
absolute latitude of the ship. 
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Figure 6.2 - Results for 12 Hour Orbiting Satellites 
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6 • 2 Time Average Population Computati ons 

For the special case of the ship at the origin of longitude and lati- 
tude, the time average population may be computed by Equation (4.3). 
Numerical integration yields a value of about 28.48% for N /n for Class II 
satellites. That is, on the average, 28.48 satellites are mutually visible 
to the ground station and the ship when they are together. 

As a comparison, it may be shown that the ratio of the area of common 
visibility Q to the area of the zonal belt A covered by the satellite 
orbits is given by 

n 

A 

when the ground station and ship are together. Hence, for Class II satel- 
lites, we obtain a value of about 22.4% for fi/A. As expected, this value 
is smaller than that for N^/N because the density f increases with latitude 
and hence contributes toward giving a higher value of in the numerical 
integration. 

The other comparison is made with the results displayed in the (0,0) 
box of Figure 6.2 which are seen to yield a smaller value than that for 

is also to be expected because the approximation made in Sec- 
tion 5.2 assumed that the satellite orbits are arcs of great circles within 
the region Q and hence yields a smaller value of the average population 
time T than that obtained by considering the actual satellite ground track. 


[l - sin (tt/ 2 - g) 
2 sin i 


( 6 . 2 ) 
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